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Abstract 

The next-to-leading open charm production in jp collisions is calculated within the 
Perturbative Fragmentation Functions formalism, to allow resummation of a s log(p^/m 2 ) 
terms. In the large px region (pr > m) the result is consistent with the fixed order NLO 
calculation, small discrepancies being found for very large pt and at the edge of phase 
space. The two approaches differ in the definition and the relative contribution of the 
direct and resolved terms, but essentially agree on their sum. The resummation is found 
to lead to a reduced sensitivity to the choice of the renormalization/factorization scale. 
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1 Introduction 



The inclusive photoproduction of heavy quarks, namely the process 

7 + H -» Q + X, 

has been evaluated up to next-to-leading order (NLO) in recent years. R.K. Ellis and P. Nason 
have first evaluated the QCD corrections to the direct jH — > QX process Their calculation 
has been successively confirmed by J. Smith and W.L. van Neerven [0]. The NLO corrections 
for parton-parton scattering, needed for the resolved part, have been produced by P. Nason, S. 
Dawson and R.K. Ellis and also by W. Beenakker et al. 0] 

More recently S. Frixione et al. have finally written computer codes || |6| capable of inte- 
grating the NLO formulas thereby producing phenomenologically useful cross sections. These 
codes will be used to perform comparisons with our results. 

Within the context of NLO calculations of heavy quark photoproduction, it is also worth 
reminding the calculation of inelastic J/ip photoproduction which has recently been presented 
by M. Kramer et al. 0. 

All the above fixed-order perturbative calculations rely on the fact that the mass of the heavy 
quark acts as a cutoff for infrared singularities. The splitting processes which involve a heavy 
quark are therefore finite order by order in perturbation theory: there is no need to subtract 
singularities when handling a heavy quark line. This has induced the authors of refs.JT], [| to 
adopt a modified version M of the MS subtraction scheme in which the heavy quark effects 
decouple in processes involving momenta much smaller than the heavy quark mass. This also 
implies that the heavy quark does not contribute to the evolution of the running coupling and of 
the structure functions, and therefore the absence of subprocesses initiated by an intrinsic heavy 
flavour generated by the structure function. All the effects of the heavy quark are therefore 
contained in the kernel cross section only. 

The NLO one particle inclusive differential distribution will however contain terms of the 
kind a s log(p^/m 2 ) which, in the large pr limit, will become large and will spoil the perturbative 
expansion of the cross section. This is reflected in a large sensitivity to the choice of the 
renormalization/factorization (r/f) scales and hence in a large uncertainty in the theoretical 
prediction. Similar potentially large logarithms, of the kind a s log(Q/m), Q being the photon 
virtuality, do appear in NLO calculations of heavy quark electroproduction total cross sections. 
The resummation of these logs is being considered in a series of papers f|. 

In Ref. [|Hj the problem of large a s log(p^/m 2 ) terms was tackled in the hadropro duct ion 



case by introducing the technique of the Perturbative Fragmentation Functions (PFF) through 
which these terms were resummed to all orders and the cross section was shown to display a 
milder scale sensitivity. 

More in detail, the PFF approach is based on the assumption that when the momenta 
involved are much larger than their mass, the heavy flavours behave as if they were massless. 
Technically, mass terms in the kernel cross sections are suppressed by powers of m/pT- This 
means that one can go back to the usual MS scheme and treat all the quarks as massless, thereby 
subtracting also the singularities of the heavy quark lines and introducing a fragmentation 



function to absorb the final state divergencies. The key point is that the massiveness of the 
heavy quark can be exploited to calculate in perturbative QCD the initial state conditions for 



the heavy quark fragmentation functions at a scale of the order of its mass 11 1. Therefore the 



PFF approach relies on the same input parameters of the perturbative calculation, but allows 
a resummation of large logarithms of (pr/m). However the neglecting of the mass terms will 
of course make the PFF result not accurate in the region px ~ Tn. In particular, the PFF 
approach is not able to calculate the total cross section, which is obtained from the fixed-order 
calculation. 



2 Scales in the photoproduction process 

As well known, when calculating the next-to-leading order correction to a direct photopro- 
duction process (i.e. the 0(ot era a 2 ) terms) one finds that singular terms associated with the 
electromagnetic vertex do appear. They are due to the emission of collinear light (massless) 
partons from the incoming photon and they are a signal for the presence of a non-perturbative 
region where the photon splits into quarks and gluons before interacting with the partons in 
the hadronic target. This problem is solved in very much the same way as the hadroproduc- 
tion case: the collinear singularity is subtracted at a given scale and factored into a structure 
function. The photoproduction process gets therefore splitted into two pieces. 

Above the aforementioned scale the photon couples with pointlike coupling to a quark which 
then undergoes a hard scattering. This term is called the point-like or direct contribution. 
Below the factorization scale the photon is regarded as a hadron, and its structure function gives 
the probability for it emitting a parton which subsequently participates in the hard scattering 
process. Hence the name of hadronic or resolved photon component. 

The photon structure functions will of course depend upon the momentum scale /z 7 at 
which the collinear singularities of the photon leg are subtracted. Neither the point-like nor 
the hadronic components are separately independent of /i 7 , because the subtracted term in the 
direct component is related to the definition of the photon structure function in the resolved 
component. 

Let us now consider the generic expression for the heavy-quark production at large px via 
fragmentation due to an on-shell photon colliding with a hadron H. Expliciting the various 
scale dependencies in the process, we write the 0(a em a 2 ) cross section in the following form 

da yH ^ Q x ~ J2 J dxdzF t H (x,/j, F )da yi ^ k (x,z, acsfan), fJ, R , fJ, F , (J,~f)D%(z, fi F ) 

ik 

+ Yl J dx i dx 2 dz F*(xi, f/ F , /j, y )F J H (x 2 , // F )d<Xy^ fc (xi, x 2 , z, a s (// K ), f/ R , fi' F , fi 7 )D^(z, p! F ) 

ijk 

(1) 

with, at the NLO here considered, 

do , 7i _ + fc(a s (// Ji ), fi R , Hf,Hj) = 



a em a s (fi R )da { ^ k + a em a 2 s (^ R )dd ( ^_ k (^ R) fx F , /i 7 ) 
a 2 s {fi R )da^ k + a 3 s (fi R )da^ k (fi R , fi F , /i 7 ) (2) 



Here fi R and /i' R are renormalization scales, /i F and /i' F are factorization scales for collinear 
singularities arising from strong interactions, and /x 7 is a factorization scale for collinear singu- 
larities arising from the electromagnetic vertex. All the structure and fragmentation functions 
obey the usual Altarelli-Parisi evolution equations 

dlog//~ Z7T Jx y 

with the exception of the photon structure function which also has an inhomogeneous 
evolution in /i 7 : 

d log 2tt % 

In the commonly-used photon density parametrizations, /i 7 is usually kept equal to so 
that the term given in eq. (§) becomes a correction to the usual Altarelli-Parisi equation (the 
so-called inhomogeneous term). 

The cross section ([!]) is of course independent of the renormalization/factorization scales 
at the perturbative order at which it is calculated. Indeed fi R cancels between a s (fi R ) and 
the explicit dependence of the kernels; fx F cancels between the structure/fragmentation 
functions and again the NLO kernel dependencies; finally, /i 7 cancels between the pointlike 
component and the photon structure function evolution in the resolved component. This last 
cancellation, in particular, prompts for the need to consider both the pointlike and the resolved 
component when evaluating photoproduction cross section. Being the factorization procedure 
entirely arbitrary, none of the two components has any physical meaning, only their sum being 
observable^]. 



x 2 + 
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0(a s 



(4) 



3 The PFF approach 



We will consider now in detail the photoproduction of heavy quarks in the framework of PFF. 
Then, following for instance ref.|l3|, we can rewrite the multidifferential cross sections for the 
photoproduction process 

1 (P 1 )+H(P 2 )^Q(P)+X 
as follows. The resolved part reads 



Ed 3 



a' 



d 3 P 



E 



dz r^-y)/* dv r 1 



ftS "777f Jl-V+VW Z 2 Jvw/z 
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1 da l \s,v) 
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5(1 - w) + s R Kij^ k (s, v, w; fj, R , n F ) 
In 



having defined the hadron-level quantities 

T 



V = 1 + 



S 



w 



-u 



S + T 



(6) 



2 For a detailed discussion on the cancellation of the ^ dependence in the photoproduction of jets, see ref. |12| 



with S = (Pi + P 2 ) 2 , T = (Pi — P) 2 e U = (P 2 - Pf. Similarly are defined the parton-level 
ones s, v e w. In terms of the momentum fractions x\, x 2 and z it holds 



XiX 2 S 



Xl 
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x 2 



zl 



(7) 



Notice that we have abandoned any distinction between the photon and the hadron factorization 
scales /i 7 and /i F . In the following we'll also usually identify the renormalization and the 
factorization scales, i.e. fi R = fi F = \i. 

The expression for the pointlike contribution amounts to rewriting (^) with the constraint 
FHx, fi F ) = 8(1 — x) with kernel cross sections for 7 scattering. By exchanging the integrations 
over v and w we can write 
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(9) 



z-VW/w' 

and we have dropped the integration limits for v, being this integration constrained by the 
^-function. 

We recall that the PFF approach to large-pr heavy quark production assumes that the heavy 
quark is produced via fragmentation of partons (charm included) which have been produced 
in a massless way in the hard scattering. The fragmentation functions can be evaluated in 
perturbative QCD at a scale of the order of the heavy quark mass and subsequently evolved 
up to the desired factorization scale \x F through the Altarelli-Parisi equations. 

The lowest order massless kernel cross sections da^^ k and dcr°^ fc (of order a 2 and a em a s 
respectively) have been long known in literature. The next-to-leading terms for massless partons 
scattering Kij_, k (of order o^) and K yi _, k (of order a em a 2 ) have been instead produced in recent 
years. Indeed Aversa et al. |14| have given the explicit expressions for the NLO massless parton- 

while Aurenche et al. [15| have instead calculated the NLO 

*k- 



parton scattering processes, Kij. 

corrections to the 7-massless parton scattering process, K 7i . 

For the numerical evaluation of eqs. (5) and (8) we have used the computer programs 
originally developed by those authors |15|] to perform the convolution of the scattering 
kernels with the structure and fragmentation functions. Some modifications had however to 
be implemented to ensure the numerical convergence of the integrals. The codes had actually 
been designed to handle light hadrons fragmentation functions. These functions typically tend 
rapidly to zero for increasing values of the argument z, and can therefore easily be numerically 
integrated. This is no more true when treating heavy quarks due to the singular behaviour of 
the Dq fragmentation function in z — 1. In ref. pi a so called "pole subtraction" procedure was 



1 



implemented. In this case we have instead modified the fragmentation function by convoluting 
it with a "regulator function" D r f(x; a, [3) given by 

D rf (x;a,(3) = A (1 - x) a x p (10) 

a and (3 have to be chosen such that D r f(x) is peaked very near x = 1 (we have taken a = 1 
and (3 = 1000), and A is the normalization factor: 

A' 1 = [ dx D rf (x;ct,0) = B(0,a + 1) (11) 
Jo 

B being the Euler beta-function. The D r f(x; a, (3) goes to zero in x — 1 fast enough to allow 
the numerical integration. We have explicitly checked that the numerical result is compatible 
- within errors - with the "exact" one given by the pole subtraction method. 

A further modification of the codes concerns the effect of the heavy quark mass m on the 
renormalization/factorization scales and on the kinematical boundaries of the integrals. Indeed 
the heavy quark mass has been taken explicitly into account by using 

flref = \jp 2 T + m 2 (12) 

as the central choice for the f/r scale and also by calculating the S, T and U invariant with 
massive kinematic. Given the heavy quark transverse momentum pt and its scattering angle 
in the center of mass frame 9 this amounts to write for the heavy quark quadrimomentum P 

t> l I Pt i 2 n PrCOsjA 

P= V~ ~a+ m 'PtA • Q (13) 
\ V sin 6 sin I 

These simple kinematical considerations have of course some clear effect in the low-pr region, 
making the numerical integration easier, but in no way allow one to reproduce the result of the 
full massive calculation, being the 0(m/p<r) terms in the kernel cross sections still missing. 



4 Numerical results and conclusions 

Let us consider the inclusive photoproduction of a c quark at present HERA energy. 

For the sake of simplicity we'll take a photon of fixed energy = 26.7 GeV which scatters 
against a E p = 820 GeV proton. This amounts to a center of mass energy of 296 GeV. According 
to the usual conventions, we'll take positive rapidities in the direction of the incoming proton. 
Then the rapidities in the laboratory frame and in the center of mass frame are connected by 
the usual formula 

1 E 

yiab = ycm+ "logTT (14) 



We use the MRSA [[TjJ structure function set for the proton and the ACFGP-mc |T7j one for 
the photon. The A^ is fixed at the central value in the MRSA set, namely 151 MeV and the 
mass of the charm quark is assumed 1.5 GeV. 

Then fig. [I] shows how the fixed order calculation (from here on referred to as FMNR) 
and the PFF one do distribute differently the overall cross section between the direct and the 



resolved contributions. The bidifferential cross section da/dp T dy is considered at the rapidity^ 
value yi ab = 1. In these plots the r/f scale is taken equal to the reference value defined in ([12]). 
Namely, \i = £/z re / with £ = 1. 

By looking at the plots we immediately appreciate that the direct and resolved contributions 
separately do behave differently in the two approaches. The FMNR direct part tend to be 
higher than the PFF one. The opposite happens in the resolved part, the PFF result being 
markedly higher, especially in the large pt region. The total cross section however, where the 
two components have been added together, displays a good similarity between the FMNR and 
the PFF results. The two curves are practically coincident up to px values around 10-15 GeV, 
and above these values the PFF one is only slightly lower. This deviation is due to the multiple 
gluon emission from the final state charm, which becomes relevant for pt 3> m. This effect is 
of course not described by the fixed order NLO calculation, which only deals with one gluon 
emission and therefore predicts an overall harder charm fragmentation function. 

The different descriptions that FMNR and PFF give of the direct and resolved part sepa- 
rately are due to the different treatment of the diagrams in fig. |3], where the photon splits into 
a cc pair. In the FMNR calculation the diagram a) is entirely included in the NLO direct part, 
the finite mass of the charm quark preventing the splitting from being singular. Consistently 
the diagram b), i.e a scattering initiated by the heavy quark, does not appear in the resolved 
part (we would have a double counting otherwise). In the PFF approach on the other hand the 
charm is massless, and the process of photon splitting has to be separated into two. Above the 
factorization scale /i 7 the contribution is assigned to the direct component, while below this 
scale the splitting process is instead considered non perturbative: it is taken into account via 
the photon structure function (fig. ||b) and assigned therefore to the resolved part. This 
different separation explains why in the PFF approach the direct component is smaller and the 
resolved one bigger than in the FMNR one. 

The direct and resolved components of the cross section are shown in fig. |2|, where the 
rapidity distribution at pt = 10 and 20 GeV is plotted. Once again, the PFF direct term 
falls below the FMNR one, due to the subtraction of part of the photon splitting to cc. By 
considering the resolved part, on the other hand, we see that FMNR displays the typical 
symmetrical bell shape usually observed in hadroproduction. The PFF shows instead a peak in 
the low rapidity region, similarly to what happens in the direct component. Still, after adding 
the two contributions the two approaches show a remarkable agreement, except at the very 
border of phase space where PFF generally produces a lower cross section. This discrepancy 
is due to the fragmentation functions being here probed exclusively at values of z near one. In 
this region they display an unphysical singularity, and resummation of Sudakov terms is needed 



to ensure a proper behaviour [11 



Next we consider the scale sensitivity of the two contributions in the two approaches. To 
this aim we recall that in the case of hadroproduction the PFF approach leads to a reduced 
sensitivity to the r/f scales, improving considerably the theoretical uncertainty at large pt- 
Fig. [| shows the uncertainty band given by varying the f/r scale between 0.25 fi re f and 2 fi re f 

3 The PFF approach actually deals with the pseudorapidity r\. In the large px limit however, where px 3> m, 
the two quantities are coincident. 



in the direct component. In both approaches £ = 2 gives the lower curve. This is coherent with 
the definition of "direct contribution" as what is above the factorization scale. Then raising 
this scale, less room is left of course for this process to take place. From the same figure we also 
notice a much larger sensitivity in the PFF approach. This is due to the fact that in this case 
a change of the factorization scale also affects the charm content of the photon. In FMNR, on 
the other hand, this process is not factorized and therefore its large contribution to the cross 
section is not affected by a variation of /i 7 . 

Fig. |5| shows the analogous band for the resolved component. For FMNR the £ = 2 choice 
still gives the lower curve. This is consistent with what was observed when studing the hadropro- 



duction case flQfl . For PFF, on the contrary, £ = 2 now gives the higher curve, at least in the 
large pt region. This is related to the fact that the photon splitting into cc is assigned to the 
resolved component below the factorization scale. 

Finally, fig. |6| shows the scale sensitivity of the overall cross section. Due to the opposite 
behaviours of the two partial contributions previously considered the two approaches can be 
seen to display similar sensitivity. 

A more detailed analysis of the scale dependence is shown in fig. [7|, where the cross section 
is plotted for three different pt values (10.5, 20.5 and 30.5 GeV) as a function of the f/r scale 
in the range £ =0.25-2. As previously noted, the dependencies are pretty similar, the PFF 
approach being slightly less sensitive. For comparison a PFF result obtained using leading 
order kernel cross sections and a s and fragmentation functions evolutions has been plotted. It 
predicts a lower cross section and displays, as expected, a larger scale sensitivity. 

To conclude, we have considered the open charm production in jp collisions to next-to- 
leading order within the Perturbative Fragmentation Functions formalism, to allow the resum- 
mation of ot s \og{p\/m 2 ) terms. Both direct and resolved components have been considered in 
detail. In the large px region (pr > m) the results are found to be consistent with the fixed 
order NLO calculations, small discrepancies being found for very large px and at the edge of 
phase space. The two approaches differ in the definition and the relative contribution of the 
direct and resolved terms, but essentially agree on their sum. The resummation is found to 
lead to a slightly reduced sensitivity to the choice of the renormalization/factorization scale. 

Note added: After this paper was written it came to our attention the work [jTSfl , where a similar 
analysis was performed. 
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Figure 1: Comparison between the fixed order (FMNR) and the Per- 
turbative Fragmentation Function (PFF) prediction of the charm pho- 
toproduction pt spectrum. 



Figure 2: Comparison between the FMNR and the PFF prediction for 
the rapidity distribution in charm photoproduction. The upper curves 
refer to pt = 10 GeV, the lower ones to px = 20 GeV. 



Figure 3: The two different ways a charm can come from a photon: via 
pointlike coupling (a) or via the photon structure function (b). 



Figure 4: Scale sensitivity of the FMNR and PFF prediction for the 
direct contribution to charm photoproduction. 



Figure 5: Scale sensitivity of the FMNR and PFF prediction for the 
resolved contribution to charm photoproduction. 



Figure 6: Scale sensitivity of the FMNR and PFF prediction for charm 
photoproduction. Both the direct and resolved component at the NLO 
are here taken into account. 



Figure 7: Scale sensitivity of the charm photoproduction cross section 
as a function of the renormalization/factorization scale. "PFF LO" has 
been obtained with leading order cross section kernels and a s and frag- 
mentation functions evolutions. 
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